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ABSTRACT 

There is not as yet full agreement on the mechanism that causes the rapid 
damping of the oscillations observed by TRACE in coronal loops. It has been 
suggested that the variation of the observed values of the damping time as func- 
tion of the corresponding observed values of the period contains information on 
the possible damping mechanism. The aim of this Letter is to show that, for reso- 
nant absorption, this is definitely not the case unless detailed a priori information 
on the individual loops is available. 

Subject headings: MHD — Sun: corona — Sun: magnetic fields — waves 
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1. Introduction 

Transverse oscillations in coronal loops have been detected in observations made with 
the EU V telescope on board of th e Tr ansition Region and Cor onal Explorer (TRACE) in 
1999 by lAschwanden et al.l (119991 ) and iNakariakov et al.l (119991 ) . Since then the detection 
of these oscillations has been confirmed and in addition damped oscillati ons have been ob 



served in hot coro nal loops by the SUMER instrument on board SOHO (IWang et al.l 12002 
Kliem et al.l 120021 ) . The TRACE oscillations have periods of the order of ~ 2 — 10 min- 
utes and comparatively short damping times of the order of ~ 3 — 20 minutes. There is 
general consensus that these oscillations are fast standing kink mode oscillations. How- 
ever, there is still debate about the mechanism that c auses the observed fast dam ping. 
Different mechani sms have been suggested: phase mixing (lOfman fc Aschwandenll2002r). res 



onant absorption (jHollweg fc Yanglll988l: 



lateral and foot-point wave leakage (jSmith et al.l Il997l : 
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to the ambient plasm a rtChen &: Schuckl 120071 ) . In order to discriminate between different 
damping mechanisms lOfman &: Aschwandenl (120021 ) suggested to study how the observed 
damping times vary a s a fu nction of the corresponding observed periods. In particular 



Of man & Aschwandenl (120021 ) claimed that the observed values are compatible with phase 



mixing if the damping time increases with period T as T 4//3 and with resonant absorption if 
it increases as T. The aim of the present Letter is to show that random observations of os- 
cillation events in coronal loops are very unlikely to produce any particular relation between 
damping times and periods whatever the actual mechanism is that causes the damping. In 
particular we focus on resonant absorption and use synthetic data for periods and damping 
times to show that various samples of pairs of periods and damping times can lead to various 
and widely different scaling laws. 



2. Analytical theory 



The suggestion that the damping time increases linearly with period, r d ~ T, was 
trigg ered by analytic a l asymptotic expressions for the damping of the quasi - mode given by 



Goossens et al.l ( 1992 ). iRuderman &: Roberts! (120021 ). iGoossens et al.l (120021 ) (see re- 



GoosscnsetaL 2006 and Goossen 
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These asymptotic expressions are derived 



e.g 
views 

in the approximation that the non-uniform layer is thin. This is the so-called thin bound- 
ary (TB) approximation. Jump conditions are used to connect the solution over the ideal 
singu l arity and to avoid so lving the non-ideal MHD wave equations (see e.g. ISakurai et al. 
199ll ; IGoossens et al.lll995l ). The asymptotic expression for the damping time in the TB 
approximation is 



Td 

T 



I c 



where ( = pij p e is the density contrast, pi and p e the constant internal and external densities, 
/ the thickness of the non- uniform layer defined in [R — l/2, R + l/2], and R the mean radius. 
F is a parameter that depe nds on the variation of density in the transi tional layer. For 



a linear variation F = A/tt (jHollweg &: Yandll988l : IGoossens et al.lll992l ). for a sinusoidal 



variation F = 2 /it ( Ruderman Sz Roberts! 20021 ) . In what follows we shall use the version 
with F = 2/tt. 

Equation (TjQ) might be considered to be an invitation for predicting a linear variation 
of damping time with period. However, the crucial point for this prediction to make sense 
is that the remaining quantities in equation (CQ) are constants; i.e. all the oscillations events 
under considerations have to occur in loops with the same density contrast ( and the same 
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radial inhomogeneity l/R. When applied to a collection of oscillating loops on which no 
prior information on their structure is available equation ([1]) does not allow any prediction. 
A naive use of equation ([T]) indicates that for a given period T the damping times can differ 
by f.e. a factor 50/3 for the combinations ( = 1.5, l/R = 0.1 and £ = 5, l/R — 0.5. 

This is indeed a naive use of equation (pQ) since T is a function of £ also. In order to make 
that clear we use results from the thin tube (TT) or long wave-length approximation for the 
period T. Our starting point is the well-known expression for the square of the frequency of 
the kink mode in a cylinder with a uniform and straight magnetic field along the z-axis (see 
e.g. 



Edwin fc Robertslll983h 



u — u k — ; > \ z ) 

Pi + Pe 

where ua = k z VA is the local Alfven frequency and va = B / ^JJ^p the local Alfven velocity 
(the subscripts "i" and "e" refer to internal and external, respectively) . Now we note that for 
the fundamental mode k z = tt/L, with L the length of the loop, and we convert frequencies 
to periods and rewrite equation (T5|) as 

— = V2A(C), (3) 

TA,i 

where ta,% = L/vA,i is the internal Alfven travel time and the function A(() is defined as 

/c + i\ 1/2 

A(C) = [^) ■ (4) 

It is convenient to adopt a reference loop with magnetic field strength B , length Lq, and 
internal density p^o- This reference loop allows us to consider loops of different dimensions 
and to introduce the normalized period T* defined as 

T 

(5) 



(V2)ta,»o 

T* can be interpreted as the period measured in the unit ta^oVZ- With this normalized 
period equation ([3]) takes the simple form 



T* = aA(Q, 



(6) 



-4- 



with a = TA,i/TA,io- From equation ([6]) we learn that, for a given value of a, T* varies 
continuously from \[2a to a when £ is allowed to vary from 1 (no loop) to oo (outside 
vacuum) . 

Let us now turn back to equation p]) and rewrite it as 



- = — B <<> 7^5- < 7 > 

TA,i VT 



with the function B(() defined as 



BiO= ( — ) 1/2 ^ = A(C)- C 1 ' 



In what follows we shall refer to the combination of the TB approximation to compute the 
damping time and the TT approximation for the period as the TTTB approximation. 



In the same way as for the period it is convenient to introduce the normalized decay 
.* 

d 



time r? as 



(2V2/7r)r Aii0 

can be interpreted as the decay time measured in the unit r^io (2v^/?r). With the use 
of this normalized decay time we can write equation (J7J) as 



rl = a B(C) -L (10) 

which tells us that, for a given value of a, depends on ( and l/R. The dependency on 
l/R is straightforward. The dependency on Q is slightly more complicated. Basically, 
is a decreasing function of ( varying from +oo for lim£ — > 1 to ^ for lim^ — > +oo. In 
summary, is a decreasing function of both l/R and 

We now combine equations (jSJ) and (jlOl . On Figure [T] we have plotted r* d versus T* for 
values of C from 1.5 up to 10 and of l/R from 0.01 to 2. For every value of a the pairs (T*, r^) 
define a vertical strip in the (T*, r^)-plane with a maximum horizontal extent of (v2 — l)a. 
As a matter of fact a can take on any positive real value, but for clarity we have limited the 
values to a = 1, 2, 3, 4, 5. The inclusion of intermediate values of a, such as a = 4/3, 5/3, 7/3, 
produces vertical strips overlapping those shown on Figure [U The iso-lines corresponding to 
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a constant value of C are vertical lines since the period is independent of the inhomogeneity 
length scale, in the thin tube approximation. The iso-lines corresponding to a constant value 
oi l/R are not straight and defined by the equation 



It is clear from Figured] that the model of resonant absorption in its most simple mathemat- 
ical formulation using the TTTB approximation produces a wide variety of combinations of 
(T*, t%). Now let us see what happens when a collection of pairs (T*, t*[) is drawn from this 
reservoir. 



The aim of this section is to show that different collections of data can be produced by 
the simple TTTB mathematical model of resonant absorption with each of them leading to 
scaling laws, r^/(T*) n = C (with C a function of the equilibrium parameters), with different 
indices, n. We use equations §6§ and (flDl) . together with the particular relation between 
t£ and T* for each index n, and derive those collections of data. As a first example we 
have plotted on Figure [2] a collection of (T*,t$) combinations drawn from the big reservoir 
depicted on Figure [TJ These combinations define a scaling law with index n = — 1. It is clear 
that we have engineered the data presented on Figure [2] so as to fit the scaling law with such 
index . Actually, the engineering is relatively straightforward. We require T*t^ = C s . The 
left-hand side is the value of the function C for a starting configuration with a = as, C = Css 
and l/R = (l/R)s- We fix a loop to start with by prescribing the values as — 1, Cs — 1-5, 
and {l/R)s = 0.1 and compute values of a, (, and l/R that satisfy the equation 



Likewise on Figure [2] we have also plotted a second collection of data that now define 
a scaling law with index n = 2. To obtain these data, we start from a loop with prescribed 
values a s = 1, (s = 5, and {l/R)s = 1 and compute values of a, C, and l/R that satisfy the 
equation 




(11) 



3. Scaling laws: how many do we want? 




(12) 




B(C) 

(MO) 



(13) 



A third collection of d ata, that now define a scaling l aw with index n = 4/3 is also shown 
in Figure [2j According to lOfman fc Aschwandenl (120021 ) this value of the index singles out 
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phase mixing as damping mechanism. Here it is obtained for a collection of data that are 
produced by the theoretical prediction for resonant absorption. We start from a loop with 
prescribed values a s — 1, Cs — 3, and (l/R)s — 0.5 and compute values of a, (, and l/R that 
satisfy the equation 



' f'\ «»\» M. m= W. (14) 



R \R) S \ a) h{CsY ^' (A(0f' 3 ' 

Finally, on Figure [2] we have also plotted a collection of data that define a scaling law 
with the canonical value n=l for the index. We start from a loop with prescribed values 
as = 1, Cs = 2.5, and {l/R)s — 0.2 and compute values of a, (, and l/R that satisfy the 
equation 



°<o. c{() = m = c+i (15) 



R \R) S c(U)' ,4«) C 

Note that a is absent from (TT5]) meaning that a solution of (fl5|) can combined with any value 
of a > 0. An obvious solution to ([TBI) is 



c = Cs, 4 = (4 1 > « > o- 



Equation ( TTdT) defines a straight line in the (T*, r^)-plane. Note that equation ( fl6l) is not the 
only solution in the (T*, r^)-plane. Any combination (T*, r^) that satisfies equation (TT5l) 
combined with a g]0, oo[ produces the same straight line in the (T*, r^) plane as solution 
( JIBjl . Since the function c(C) is a decreasing function of its argument, a value C > Cs is 
combined with a value l/R < (l/R)s and conversely a value C < Cs is combined with a value 
Z/i2 > (l/R) s . 

It is clear from Figure [2] that almost any scaling law can be obtained from data produced 
by the simple TTTB mathematical model of resonant absorption. All the periods and 
damping times plotted in Figure [2] correspond to combinations of C and I / R, obtained from 
equations ( fT2i) . ( Tl3j) . (tHJ), and ( |T5l) . that are reasonable values of these quantities in the 
ranges C € [1-5, 10] and 1/ Re [0.016, 1.59]. The discrete sets of values in Figure [2J arise due 
to the considered discrete values of a, but loops with different internal travel times, with 
respect to the reference loop, would give the full set of values for and T* along the solid 
lines depicted in Figure |2j 
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4. Beyond the TTTB approximation 

The TTTB approximation as a mathematical model for resonant absorption has its clear 
limitations. First, the values of the period defined by equation (jfJJ) are independent of the 
radius R and the inhomogeneity length scale l/R. In reality, the period is a function of the 
density contrast ( , l/R, and R. As a consequence the iso-lines corresponding to a constant 
value of ( in the (T*, r^)-plane are no longer straight vertical lines. For a given value of a 
this dependency of the period on ( produces additional spread on the original vertical strips. 
Second, the TTTB approximation is an accurate approximation as long as the damping time 
is sufficiently longer than the period, since this is an inherent assumption for carrying out the 
asymptotic analysis leading to equation flTJ. We can expect this equation to be inaccurate 
for large values of the density contrast ( and large values of the inhomogeneity length scale 
I 1 R. Figure [3] is the twin versio n of Figure [J with the values of T* and now computed 



for fully non- uniform 1-D loops (IVan Doorsselaere et al.l 12004 ). Figure [3] is very similar in 



appearance to Figure [TJ The differences occur where they expected to occur. The vertical 
strips of Figure [T] are now replaced with strips slightly oblique to the vertical axis and the 
largest differences appear at the short values of the damping times corresponding to the large 
values of the density contrast ( and large values of the inhomogeneity length scale l/R. The 
basic message from both figures is the same. The model of resonant absorption produces a 
wide variety of combinations in the (T*, r^)-plane. If we draw a curve defined by a given 
relation between T* and t$, for instance r^/(T*) 2 = C, in this plane we can graphically 
determine as many points on this curve as we want. If we select the data determined in this 
way the result is a scaling law with index n = 2 as the one in Figure [2j The only difference 
is that now the procedure is numerical all the way. 



5. Conclusion 

In this Letter we have explained that scaling laws for a group of observations of oscil- 
lating loops to discriminate between different damping mechanisms are of not much use if 
there is no a priori information on the properties of the loops. The analytical expressions 
obtained by the TTTB approximation enable us to show in a straightforward and easy to 
follow manner that in the framework of resonant absorption collections of synthetic data can 
be produced that follow almost any scaling law. Then, we have backed up our findings by 
numerical eigenvalue computations that do not suffer from the TTTB limitations but are on 
their own rather less instructive. 

It might be argued that nature is not as cunning as the authors of this paper and does not 
attempt to engineer data according to prescribed rules as those defined in these equations. 
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On the other hand, there is no reason why nature would want to produce coronal loops that 
all have the same values of ( and l/R. So, unless there is accurate a priori information on 
the coronal loops available, the use of scaling laws to discriminate between different damping 
mechanisms is questionable, to say the least. 
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Fig. 1. — Normalized damping time as a function of the normalized period for discrete values 
of C e [1-5, 10] and l/R G [0.01, 2], and five values of a. 
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Fig. 2. — Four synthetic scaling laws corresponding to n — —1,4/3,1,2, and values of a, 
from left to right, a = 1, 2, 3, 4. 




Fig. 3. — Normalized damping time as a function of the normalized period for fully non- 
uniform loop models with ( G [1.5, 20] and l/R e [0.1, 2]. The six vertical strips correspond 
to six values of k z = n/L. 



